Optically trapped quasi-two-dimensional Bose gases in random environment: quantum 

fluctuations and superfluid density 
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We investigate a dilute Bose gas confined in a tight one-dimensional (ID) optical lattice plus a 
superimposed random potential at zero temperature. Accordingly, the ground state energy, quantum 
depletion and superfluid density are calculated. The presence of the lattice introduces a crossover to 
the quasi-2D regime, where we analyze asymptotically the 2D behavior of the system, particularly 
the effects of disorder. We thereby offer an analytical expression for the ground state energy of 
a purely 2D Bose gas in a random potential. The obtained disorder-induced normal fluid density 
n n and quantum depletion n<j both exhibit a characteristic 1/ In (1/7220020) dependence. Their 
ratio Un/rid increases to 2 compared to the familiar 4/3 in lattice-free 3D geometry, signifying a 
more pronounced contrast between superfluidity and Bose-Einstein condensation in low dimensions. 
Conditions for possible experimental realization of our scenario are also proposed. 
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I. INTRODUCTION 

The effect of dimensionality of a bosonic system on 
the presence and nature of the Bose-Einstein condensa- 
tion (BEC) as well as on the superfluid phase transi- 
tion has received long-standing interests both experimen- 
tally and theoretically 0, . The physics at low dimen- 
sions exhibits fundamental difference from that in three- 
dimension (3D). In particular, the strong long-range 
phase fluctuations typical of low-dimensional bosonic sys- 
tems usually inhibit the formation of long-range order, 
which on the other hand characterizes the 3D BEC and 
corresponding phase transition at low temperature 0]. 

Earlier work on low-dimensional bosonic systems 0] 
have culminated in, particularly in uniform 2D case, two 
important theoretical discoveries. The first is that in 
2D a true condensate can only occur at T = and its 
absence at finite temperature follows from Bogoliubov 
/c _2 _J4[ or Hohenberg-Mermin- Wagner (BHMW) theorem 
0> HI- On the other hand, a superfluid phase transition 
has been proven to exist at sufficiently low temperature in 
2D 0,11]. However, according to Kosterlitz and Thouless 
(KT) [9( , such transition is associated with the unbinding 
of vortex pairs or quasi-long-range order, in contrast to 
the 3D phase transition that features long-range order 
parameter. Below the KT transition temperature, a 2D 
Bose gas (liquid) is characterized by the presence of a 
"quasicondensate" [l(| HU • 

The remarkable experimental progress with ultracold 
atomic gases, especially in cooling and confining of cold 
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atomic gases in traps with controllable geometry and 
dimension, have significantly stimulated new interests 
in low-dimensional systems [12j . Tight confinement in 
one or two directions considerably affects the properties 
of Bose gases such as collisions and phase fluctuations 
fl3| , introducing a crossover to the quasi- low-dimensional 
regime. As such, quasi-2D quantum degenerate Bose 
gases have been experimentally produced both in single 
"pancake" traps and at the nodes of ID optical lattice 

However, these marginal 2D Bose gases are qualita- 
tively different from corresponding infinite ones. Along 
this line, Petrov et al. [l3[ have pointed out that the pres- 
ence of trapping potential suppresses long-range thermal 
fluctuations and that in a quasi-2D system a true conden- 
sate can exist within a wide parameter range. This the- 
oretical prediction has been echoed by Fischer et al. [IU 
who obtained in a marginal 2D case a model-independent 
geometrical equivalence of the BHMW theorem. 

Compared with harmonically trapped systems, opti- 
cally trapped Bose gases allow more experimental con- 
trollability with tunable inter-atomic interactions, tun- 
neling amplitudes between adjacent sites, atom filling 
fractions and lattice dimensionality [l2| , thereby present- 
ing a more useful testing ground for theoretical ideas 
in studying low-dimensional systems in novel conditions. 
On the other hand, disorder has been observed to cause 
dramatic influence on a BEC and has attracted huge in- 
terests recently [H, [l?} • In view of the availability to con- 
trol a ID optical lattice and external randomness, there- 
fore, one especially appealing direction of investigation 
consists in studying the effect of external randomness on 
a Bose gas trapped in a ID optical lattice. 

In this paper, we investigate the ground state proper- 



2 



ties and superfluidity of a lD-optical-lattice trapped Bosc 
gas in a random environment at T = 0. Capitalizing 
on the characteristic lattice-induced 3D to quasi-2D di- 
mensional crossover, we analyze effects of disorder in the 
asymptotic 2D regime. The present work is composed 
of two parts. In the first part, we calculate the ground 
state energy and quantum depletion for the model system 
using the path integral approach within the Bogoliubov 
approximation. Discussion on the dimensional crossover 
property in a random potential is presented. In partic- 
ular, our results in the quasi-2D regime with varnishing 
disorder are in good agreement with that of a homoge- 
neous 2D Bose gas at T = (a Il8j|. We suggest, there- 
fore, that our result gives the analytical expression for 
the ground state energy of a uniform 2D dilute Bose gas 
in the presence of weak disorder. In the second part, we 
calculate the disorder and lattice induced normal fluid 
density n n at T — 0. Our results in the anisotropic 
3D regime reproduce the well known ratio n n jn^ = 4/3 
[l~9l |20| with rid being the quantum depletion due to dis- 
order. Whereas, in the quasi-2D regime, n n exhibits a 
1/ln (l / V^d^Id) dependence unique of a 2D system and 
the ratio becomes asymptotically n n /rid = 2, indicating 
a more pronounced contrast between superfluidity and 
BEC in low dimensions. 

The outline of the paper is as follows. In Sec. II, 
we introduce the grand canonical partition function for 
a dilute Bose gas in the presence of a ID optical lattice 
and weak disorder at T = 0. Accordingly, the analyti- 
cal expressions for the ground state energy and quantum 
depletion are derived. Sec. Ill presents a detailed dis- 
cussion on the dimensional crossover in the ground state 
properties induced by a ID optical lattice. Effects of dis- 
order in the crossover regimes are analyzed. In Sec. IV, 
we calculate the supcrfluid density and study its behav- 
ior respectively in the 3D and quasi-2D regime. Finally, 
we summarize our results in Sec. V and propose possible 
experimental scenarios. 



II. BOSE GASES IN THE PRESENCE OF A ID 
OPTICAL LATTICE AND WEAK DISORDER 

A. Path integral approach 

Our starting point is the grand-canonical partition 
function of a 3D weakly interacting dilute Bose gas [2l[ 
in the presence of a ID optical lattice and weak disorder 



£>[V*,V]e ~- 



where the action functional S [ip* , ip] reads 



(1) 



S[i>*,i>] = J dr I dr^*(r,r) 



a a 2 v 2 

dr 2m 



+ V opt (r) + V ran (r) + ^(r,T) 



In Eqs. ([J) and [-0* (r, r) , tp (r, r)] collectively de- 
note the complex functions of space and imaginary time 
r, ft = 1/fcgT with ks being the Boltzmann constant 
and T being the temperature, /i is the chemical potential 
and g e is the effective two-body coupling constant in the 
presence of a ID optical lattice. The V op t(r) and V ran (r) 
respectively represent the ID optical lattice and external 
random potential. 

The optical potential V op t (r) in Eq. © is given by 



V opt {v) = s x E R sm 2 (q B z), 



(3) 



where s is a dimensionless factor labeled by the intensity 
of laser beam and En = T?q 2 B j2m is the recoil energy 
with hqs being the Bragg momentum. The lattice period 
is fixed by qs = 7r/d with d being the lattice spacing. 
Atoms are unconfined in the x — y plane. 

Disorder V ran (r) in Eq. ^ is produced by the random 
potential associated with quenched impurities [l9l. l20H22l] 



V dts (v) 



N imp 

E 

i=i 



u(|r-ri|), 



(4) 



V(r,r). (2) 



with v(r) describing the two-body interaction between 
bosons and impurities, being the randomly distributed 
positions of impurities and Ni mp counting the number of 
rj. Here, we restrict ourselves to the conditions of a dilute 
BEC system in the presence of a very small concentration 
of disorder. Thereby, v(r) can be approximated by an 
effective pseudo-potential in the form v(r) — gi mp 5(j) 
[l9l | . with Qi mp being the effective coupling constant of 
an impurity-boson pair confined in a ID optical lattice. 

It's important to mention that the tight confinement 
in the direction of optical lattice considerably influences 
the value of effective coupling constant [13|, [23[ in Eq. 
([2]). Particularly, in the presence of optical lattice, g e 
generally exhibits dependence on the density and lattice 
parameter [24j , in marked contrast to a free 3D Bose gas 
where g^o = 47rfi 2 a3£)/m with 0313 being the 3D scatter- 
ing length. For formulation clarity, however, below we 
shall use g e and gi mp for notational convenience while 
leaving aside their specific expressions in order to ob- 
tain general expressions for the ground state energy and 
quantum depletion. Analysis of the lattice-renormalized 
effective coupling constant will be given in Section IV. 



B. Beyond- mean-field ground state energy and 
quantum depletion 

In what follows, we focus on the situation where the 
optical lattice is strong enough to create many separated 
wells that give rise to an array of condensates; while full 
coherence is still assured by the quantum tunneling. By 
this assumption, one can refer to hq as the condensate 
density and neglect the Mott insulator phase transition. 
We also suppose disorder is sufficiently weak. Under 
these conditions, one is able to investigate the ground 
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state properties of the model system using Bogoliubov's 
theory (jf. 

We shall restrict ourselves to the case where s is 
relatively large that the interwell barriers are signifi- 
cantly higher than the chemical potential [i [25j |. We 
thereby only consider the lowest Bloch band where the 
condensate, in the tight-binding approximation, can be 
written in terms of Wannier functions as <pk z (z) = 
Y,ie Uk: "w(z - Id) where w{z) = exp(-z 2 /2a 2 )/n 1 / 4 a 1 ^ 2 
with d/a ~ 7TS 1 / 4 exp(— 1/4-y/i). Expanding the bosonic 
field variables in Eq. ([5} by the expression ip (r, r) = 
Ek,n^,nfe(z)e- l(tx+i » 9) e !U " T with w„ = 2nn/hf3 be- 
ing the bosonic Matsubara frequencies where n are inte- 
gers, the action Eq. ^ takes the form 



npv 



k , n 



+ ~2 V'k+q,n+mV'k'-q,n'-m , /'k',n'V'k,' 

k,k',q 

n,n ,m 

+ X! ^k-k'V'k.nV'k'.n- 



(5) 



k,k' ,n 



Here, 



(h 2 /2m)(kl + kl) - 2i[l - cos(fc a d)], with i 



being the tunneling rate between neighboring wells, is the 
energy dispersion of the noninteracting model, V is the 
volume of the system and g e is the lattice renormalized 
coupling constant given by 




x =2t/g e n 



x = 2t /Jr e "o 



FIG. 1. (a): Scaling function F(x) in Eq. ([10)1 (solid line) and 
its asymptotic behavior (dashed line), (b): Scaling function 
H(x) in Eq. (|13|l (solid line) and its asymptotic behavior 
(dashed line). 



we set u = ffe^o to ensure a gapless quasiparticle spec- 
trum [27| . 

In the continuum limit, the sum in Eq. ([7]) is re- 
placed with integrals. To avoid the large- k divergence 
in the integration over k x and k y , however, one must in- 
troduce a renormalization of the coupling constant by 
replacing g e -> g e - (g 2 e /V) J2k^o ( 1 / 2e k) and 5™ P ~> 

~9*m P - (gLp/V) Ek/p (1/24) in Eq. ©. Upon this 
replacement, one obtains after integration 



9t 



d/2 

d I w (z)dz 

d/2 



9t 



2ira 



(6) 



In Eq. ([5]) , the Vfc is the Fourier transform of V ran (r) = 
Y^i9impS( r ~ r d witn 9im P = g imp d/ \f2~na being the 
lattice-renormalized impurity-boson coupling constant, 
i.e. Vk = (1/V) J e 4k r V r ran (r)dr. For simplicity, the ex- 
ternal randomness is assumed to be uniformly distributed 



with density 



Nimp/V and Gaussian correlated 



[26l | . Hence the two basic statistical properties of the 
disorder are the average value (Vq) = gi m pni mp and the 



correlation function (V^kVk) = g 



,-,2 



imp n irnp 

notation (..) stands for the ensemble average over all pos- 
sible realization of disorder configurations. 

By applying the Bogoliubov theory to the action ([5]) 
and proceeding in the standard fashion [3j, one obtains 
the zero-temperature thermodynamic function f2 = E g — 
V[ino with the ground state energy E g reading 



/V. Here, the 
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i 



~9eTi\ 



2V?' 

k#0 



g e n 



n 



~2 

Wimp!]; 



V 



imp x ^ 

^E 2 

k^O k 



(7) 



(7g7iQ is the energy 



Here, £ k = V( e k - M + 2 <?e?V 
spectrum of the elementary excitations and no is the con- 
densate density. In conformity with the general theory, 



mRg e 
2irh 2 d 



arccoth 



2t 

g e n 



>(8) 



where the two parameters 7 = 2Kgi mp /g e with k 
rii m p/no and 



R 



^imp ^9imp 

m g 2 



(9) 



characterize the strength of disorder in a ID optical lat- 
tice. In Eq. (jHJ, the function F(x) with the variable 
x = 2tj (g e no) is defined as 



F(x) 



(x + 1) 



(3x + 1) arctan 



— 3\/x 



-I In 
2 



2a; + 1 + 2^2: (x + 1) 

tan _1 (2) 



7rarcsinh (v^) + 2 / 
./o 



<iz. 



(10) 



The integration in Eq. (| 10[) can be easily performed 
numerically and the result is shown in Fig. [IJa). In the 
ground state energy Eq. ((SJ) , the first two terms give the 
mean-field contribution modified by an optical lattice and 
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disorder; whereas, the last two terms represent beyond- 
mean-field corrections, as a consequence of quantum fluc- 
tuations respectively induced by interatomic interaction 
and external randomness. 

Quantum depletion (AN = N — No) refers to the aver- 
age number of atoms with nonzero momentum Q which 
can be calculated within the Bogliubov's theory as 



k#0 



e k + 9en -E k 
2E U 



~2 V°k 
n n imp9imp 



.(ii) 



By replacing the sum with the integral in the continuum 
limit, one obtains 



AN 



mg e 



2it 2 h 2 d 



H 



2t 

g e n 



ttR 
~8~ 



2t 

g e no 



(12) 

where the function H(x) with x — 2tj (g e no) is defined 
as 



H(x) = (x + 1) arctan ( — — 



(13) 



III. 



DIMENSIONAL CROSSOVER FROM 3D TO 
QUASI-2D AND 2D REGIMES 



At low energies, the physical properties of a dilute Bosc 
gas can be expressed in terms of the two-body scattering 
amplitude [2g|. It has been well established that a tight 
confinement along one or two directions will considerably 
affect the scattering properties of atoms, particularly, in- 
troducing a dimensional crossover from anisotropic 3D to 
low-dimensional regimes [13l . [24j . 

The two-body scattering problem in the presence of 
a ID optical lattice has been analytically investigated 
in Ref. [Hj]- For sufficiently deep lattices and chemical 
potential fi, which is small compared to the interband 
gap, two distinct regimes can be identified: (i) for /i <C At, 
where the wavefunction spreads over many lattice sites, 
the system retains an anisotropic 3D behavior. In this 
limit, Eq. takes the limiting form g e = g^o with 



33D 



Anh z a iD 



(14) 



with a^jj = a^pd/ (\J2tt(j) being the lattice-renormalized 
s-wave scattering length; (ii) for \i ^> At, the tunnel- 
ing between wells is negligible, and the two interacting 
bosons are in the ground state of an effective harmonic 
potential characterized by frequency ujq = h/ma 2 and 
harmonic oscillator length a. In this limit, the system 
undergoes a crossover to the quasi-2D regime where the 
coupling constant is reduced to that in a tight confined 
harmonic trap g e = gnd [HI, [U 29j where 



2J2^h 2 



1 



gh 



Til 



0-2D/o,ZD 



(1/V2tt) In [l/n 2D a 2 D ] 



(15) 



with the surface density n 2 £, — npd and the effective 
2D scattering length a 2D = ^h/mujQ = a [13l |. With 
decreasing a, the 2D features in the scattering of two 
atoms become pronounced [l3j . And in the limit a <C a, 
Eq. (|15p becomes independent of the value of 0315 and a 
regime of purely 2D scattering is achieved with Eq. (fT5|) 
reducing to the coupling constant of a purely 2D Bosc 
gas g h -> g 2 D where 



.9213 



AttK 2 



m \n(l/n 2 Da\ D ) 



(16) 



Here the logarithmic dependence on the gas parameter 
n 2 Da 2 D i s unique of the 2D geometry. 

Taking into account of the dimensional crossover in 
the effective coupling constant, below we focusing on an- 
alyzing the behavior of the ground state energy in Eq. 
© and quantum depletion in Eq. (IT21 . respectively 
in the anisotropic 3D and 2D geometry. In the limit 
2t/nog e 1, corresponding to the anisotropic 3D regime, 
we find F(x) ~ 32/15-y^, as is shown in Fig. QJa) with 
the dashed curve. Substitutions of this limiting value in 
Eq. ([8]) together with Eq. (IT4l) . yield the ground state 
energy of an effectively free 3D Bose gas composed of 
bosons with effective mass m* = h 2 /2td 2 and coupling 
constant g^o [23] 



E g l_ 2 

T = 2 93Dn ° 



~T \ 1/2 

n a* D 




(17) 



In Eq. (fTTj) . the two characteristic parameters of dis- 
order in Eq. ([S]) respectively take their 3D value, i.e. 
7 = K&3£>/a3D an d and R3D = K i>^ D /a\ D , showing the 
3D feature of the interaction between the impurity-boson 
pair. The first term in Eq. (1171) represents the mean- 
field ground state energy; whereas the remaining terms 
exhibit the f amiliar dependence on the effective 3D gas 
parameter nod^ n , thereby consisting of the generalized 
LHY correction [32| to the presence of a ID optical lattice 
and weak disorder. Eq. (|17[) bears formal resemblance 
with the corresponding result in Ref. [2(| which deals 
with a 2D optical lattice system, in consistent with the 
effective mass theory in the 3D limit where the lattice 
system is effectively treated as a free gas with effective 
mass and coupling constant. The main difference is re- 
lated to the value of renormalized coupling constant g^u 
where the renormalization factor for different lattice di- 
mensions [30| . 

In the opposite 2D regime where 2t/g e riQ 1 and a <C 
a, F(x) exactly approaches a limit F(x) — ■n/A—n/2 log a; 
with logx ~ \n(mt / n 2 r>2Trh 2 ) + In [in (l/^^a;^)] , as 
shown in Fig. [Ua) w ith the dashed line. In this limit, 
the Bloch dispersion can be neglected and the scattering 
problem reduces to 2D with the coupling constant Eq. 
(ITBl . In such conditions, Eq. ([5]) yields the ground state 
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energy of a 2D Bose gas in the presence of disorder 



E, 



g2D 



L 2 



1 2 
-j9-2DU w 



In [in (l/n 2 Ddl D )] 
In (l/n 2 z?a| D ) 

arccoth 



72D + 2i?2L>- 



1 



2t 



ln(l/n 2D a 2 D ) 



where L 2 is the surface area of gas, n 2 D = n^d is the 
surface density and B = 1/2 — In (mt/n 2 D2nh 2 ) . In ad- 
dition, the two parameters of disorder respectively take 
their 2D value 72D and R 2 d- Both parameters, how- 
ever, depend the 2D expression of ft mp which needs to 
be obtained from investigating in detail the 2D scatter- 
ing problem of a boson with a quenched impurity. Such 
problem is definitely non-trivial, and shall be left for fur- 
ther investigation in the future. In spite of this, Eq. (fT5|) 
has shed light on the ground state properties of a 2D 
Bose gas in the presence of weak disorder. 

Particularly, Eq. (JTBJ) presents one of the key results 
of this paper as follows: First, Eq. (fTS)) in the absence 
of disorder formally reproduces corresponding results in 
Ref. [29j for the ground state energy of a purely 2D dilute 
Bose gas. From this viewpoint, we expect that the char- 
acter of a lD-lattice-confined Bose gas in the presence 
of weak disorder in the 2D regime will be similar to a 
purely 2D Bose gas in a random potential. Therefore, we 
argue that Eq. (fT5|) provides an analytical expression for 
the ground state energy of a uniform 2D Bose gas in the 
presence of weak disorder. Specifically, the last two terms 
provide the contribution of disorder to the ground state 
energy. Second, Eq. (TTg|) has provided beyond mean-field 
corrections due to quantum fluctuations in the 2D geom- 
etry. These corrections arise from combined effects of 
interatomic interaction and disorder, and exhibit in 2D a 
characteristic 1/ In (l/n 2 Dd 2D ) dependence, in contrast 

to the 3D counterpart \J noa'^ D . 

In a similar fashion, we analyze the asymptotic behav- 
ior of quantum depletion. In the limit 2t/g e no 1, 
corresponding to the anisotropic 3D regime, H(x) ~ 
2/(3y/x). Consequently, one finds the quantum deple- 
tion in 3D, 



good agreement with Ref. [lg, [29| on the quantum de- 
pletion of a purely weakly interacting 2D Bose gas. The 
B second term in Eq. (|2U|) . therefore, presents the disorder- 
In (^/n 2 Da 2D ) induced condensate depletion in 2D. Furthermore, com- 
parison of Eq. (f20f with Eq. (fT9)l shows that, in the 
region where Bogoliubov theory applies, for the same 
(18) value of the gas parameter the quantum depletion due 
to disorder is larger in 2D than in 3D. Similar conclusion 
has been drawn in Ref. [29j for the quantum depletion 
induced by interatomic interaction. 



AN 



o 7r - \ rrr 



-o . 1/2 

n a% D x 



(19) 



characterized by the dependence on the 3D gas parameter 
(iioa^jj) 1 / 2 . In the opposite 2D limit, on the other hand, 
9e = 92D<i and H(x) saturates to the value 7r/2. Eq. 
(fT2| thereby asymptotically approaches the 2D quantum 
depletion as 



AiV 



2D 



R 



2D 



1 



In (l/n 2D a 2 D ) '' 



(20) 



IV. SUPERFLUID DENSITY 

In this section, we calculate the superfluid density of a 
dilute Bose gas in the presence of a ID optical lattice and 
weak disorder. The general definition of the superfluid 
density is proposed by Hohenberg and Martin [f|. We 
emphasize that superfluidity is a kinetic property of a 
system and superfluid density is essentially a transport 
coefficient, in contrast to the condensate density which 
is an equilibrium quantity. Superfluid density can be 
determined by the response of the system to an external 
perturbation |6j. 

In this paper, we adopt following definition: supposing 
that a linear phase Q • r is imposed on the originally static 
bosonic field which gives rise to a superfluid velocity v — 
hQ/m; in response, the thermodynamic potential of the 
system is changed by [33l435j 



sn _ h 2 ^ 
V 2m ^— ' 

ap 



(21) 



Here, the transport coefficient n ap - is interpreted as the 
superfluid density [35j |. In general, the n a p is a tensor for 
an anisotropic system. 

To obtain n Q ^, we substitute the wavefunction for a 
flowing condensate i/j (r, r) = (p (r, r) e lQ r into Eq. ( [2]) 
and obtain the action Sq for the superfluid 



k.n 



h 2 

/kQ + 

2m 



^k,n, (22) 



where S refers to the action for a static BEC in Eq. 
©, and / k q = [h 2 /m (k x Q x + k y Q y ) + 2Q z tdsin (k z d)] . 
Proceeding in a similar fashion as in Sec. II, we obtain 



9enl 



£!q = V(-p,no 



k^O 

~2 \ " 

n imp9imp n / 

k^O 



£k - A* + ft™0 



El 



ft 



(23) 



kQ 



which is proportional to 1/ In (I/^dc^d) > the small pa- 
rameter in 2D. For varnishing disorder, Eq. (I20[) is in 



where 
though fx - 



= Vfck M + 2 5e«o) 

Li-H 2 Q 2 /2m. 



ni depends on Q 
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Since the presence of a ID optical lattice breaks the 
global rotational symmetry and leaves the gas system 
only isotropic in the x — y plane, one can write n a /3 = 
n aa 5 a j3 where n xx = n yy ^ n zz . Expanding Eq. ([25)1 
in powers of Q and truncating at the quadratic order, 
we compare the resulting expression with Eq. (|21j) and 
obtain 



2n imp g 2 n ^ h 2 k 2 x e£ 



V 



k#0 



m 



and 



2mni mp gf no ^ e£ _ s 
" Z-^i [2tdsm(k z d)] 

k#0 k 



h 2 v 



(24) 



.(25) 



Similar results have been obtained in Ref. [20j using 
current-current response function. The formal agreement 
between the two affirms that, in spite of different ways 
to impose perturbation and various options of physical 
quantities to measure the response, these different routes 
to obtain superfluid density can be unified within the 
framework of the linear response theory. 

The disorder-induced normal fluid density fraction can 
be obtained through (n n ) a g = (1 — n aa /n) 5 a p. Taking 
the continuum limit of Eqs. (|24l) and (1231) . one finds 



ij^n)xx iV'Vjyy R 



5 m 9e 



2t 



8h 2 ird \g e no J ' 



(26) 



and 



(n n )zz = R 



/ m \ 2 1 
Vm* / 167rn 



K 



m* ) 167rnod 3 \g e no 



2f 



(27) 



where I(x) and K{x) are functions of variable x 
2t/g e riQ respectively defined as 



/(*) 



vT 



— x In 



1 + VT 



(28) 



and 



K [x) = In 



1 + VT 



(2 - jb) VT+~: 



,(29) 



The results of Eqs. J28J) and ([29]) are plotted in Fig. 
[2j In the asymptotic 3D limit, one finds I(x) ~ 2/3y/x 
and K{x) ~ 4/3a: 3 / 2 , corresponding the dashed curves in 
Fig. [21 In such situation, Eqs. (|2^|) and ( 1270 respectively 
become 



2lT ~ , ,,, 




and 



2ir ~ I m 
{n n )zz — —R3D\ — - 
3 V m 



n al D 



(30) 



(31) 



Eq. (|3"Tj) demonstrates similar dependence on the 3D gas 
parameter as the 3D quantum depletion in Eq. (1191) . 




x=2t/g e n 



FIG. 2. (a): Scaling function I(x) in Eq. ([28]) (solid line) and 
its asymptotic behavior (dashed line), (b): Scaling function 
K(x) in Eq. (1290 (solid line) and its asymptotic behavior 
(dashed line). 



Moreover, the ratio between Eq. ([3TJ)) and the disorder 
induced quantum depletion rid in Eq. (fl9)) equals 4/3 
in the unconfined X(Y)direction, in agreement with Ref. 
(l9| ; whereas, this ratio becomes (n n ) zz /nd — im*/3m 
due to the increased inertia of the gas along the direction 
of optical lattice [2(| ■ 

In the opposite 2D limit, one obtains the limiting ex- 
pression I(x) ~ 1 and (n n ) 2 D = (n„) x; 
to be 



(n n ) yy is found 



K) 



R.2D 



1 



2D 



2 ln{l/n 2D al D ) 



(32) 



Equation ([3"2"j) presents another key result of this paper, 
providing an analytical expression for the normal fluid 
density in a homogenous Bose fluid in 2D in the presence 
of weak disorder. Eq. (13"2l shows that the normal fluid 
density in 2D exhibits a characteristic 1/ In (l/^zja^) 
dependence. With respect to the 3D case, a comparison 
of Eqs. (|20[) and (|32|) leads to n n /rid = 2 in 2D, indicating 
a more pronounced contrast between superfluidity and 
Bose-Einstein condensation at T = 0. On the other hand, 
K(x) in Eq. (l2T)t diverges in the limit x — > 0, leading to 
diverging n zz Eq. (|270 for vanishing tunneling. This 
signals the absence of superfluidity along the direction of 
optical lattice, which is consistent with the kinematical 
2D nature of the Bose gas in the absence of tunneling 
along the direction of the laser. 



V. POSSIBLE EXPERIMENTAL SCENARIOS 
AND CONCLUSION 

Central to testing the validity of the physics in this 
article concerns experimental realization of a BEC in 
the superfluid phase along the entire evolution from 3D 
to quasi- 2D. Present facilities have allowed one to ad- 
just the depth of an lattice, realize tight confinement of 
the motion of trapped particles and ultimately achieve 
a kinematically 2D gas. In typical experiments to date, 
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quasi-2D quantum degenerate Bose gases have been ex- 
perimentally produced both in single pancake traps and 
at the nodes of ID optical lattice potentials [l4[. In ad- 
dition, it has been suggested that BEC and superflui dity 
can be both achieved below a critical temperature [13j . 
Furthermore, adding a tunable periodic potential allows 
one to combine the benefit of the reduced dimensionality 
with the advantage of working with large yet coherent 
samples [25j . 

Upon overcoming above difficulties, the experimental 
realization of our scenario amounts to controlling three 
parameters whose interplay underlies the physics of this 
work: the strength of an optical lattice s, the interaction 
between bosonic atoms grig , and the strength of disorder 
R. All these quantities are experimentally controllable 
using state-of-the-art technologies. The interatomic in- 
teraction can be controlled in a very versatile manner 
via the technology of Feshbach resonances (36[. In the 
typical experiments to date, the values of ratio grig/En 
range from 0.02 to 1 [37U38|. The depth of an optical lat- 
tice s can be changed from OEr to 32Er almost at will 
39]. Disorder may be created in a repeatable way by 



introducing impurities in the sample 40 1. or using laser 
speckles and multi-chromatic lattices 



41 



Further difficulties may arise in measuring the beyond- 
mean-field corrections to the ground state energy along 
the dimensional crossover. For typical values of the atom 
density and scattering length, such corrections remain 
very small and hard to observe in usual experiments that 
measure density profiles or release energy. They can be 
visible, however, in the frequencies of collective excita- 
tions in a lattice system [2a, 0, 5j| . The direct mea- 
surement of quantum depletions of a quasi-2D conden- 
sate can be achieved either through observing ballistic 
expansion or applying Bragg spectroscopy (47). It 
is worth mentioning that the possibility to use ballistic 
expansion to measure quantum fluctuations is associated 
with the characteristics of an optical lattice where the 



confinement frequency at each lattice site far exceeds the 
interaction energy. As such, the time-of-flight images are 
essentially a snapshot of the frozen-in momentum dis- 
tribution of the wavefunction at the time of the lattice 
switch-off, thus allowing for a direct observation of quan- 
tum depletions. This technology cannot be applied, for 
example, to measure quantum depletions of a quasi-2D 
Bose gas confined in a harmonic trap. From this per- 
spective, Bragg spectroscopy admits broader ranges of 
application, independent of methods of confinement to 
create quasi-2D BEC's systems. 

We expect, therefore, that the phenomena discussed 
in this article should be observable within current exper- 
imental capability. We emphasize here that the presented 
work is restricted to weak disorder and weak interatomic 
interaction. For further investigations in the presence 
of stronger interatomic interaction or disorder, the path- 
integral Monte Carlo simulation is a reliable method [48[ . 

In summary, we have investigated a dilute Bose gas 
trapped in a ID optical lattice and a random potential. 
Capitalizing on the characteristic dimensional crossover 
properties, the obtained results in the quasi-2D regime 
allow us to derive analytical expressions for the ground 
state energy, quantum depletion and superfluid density 
of an effectively pure 2D Bose gas in the presence of weak 
disorder. Our analysis signifies a more pronounced effect 
of disorder in systems with reduced dimensionality in en- 
hancing quantum fluctuations and depleting superfluid 
density. In particular, the ratio between the normal fluid 
density and the corresponding condensate depletion in- 
creases to 2 in 2D, in contrast to the familiar 4/3 in 
lattice-free 3D geometry. 
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